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ABSTRACT 

We present the results from an ultra-high-resolution 7 mm Very Long Baseline Array (VLBA) study of the 
relativistic jet in the BL Lacertae object OJ287 from 1995 to 2011 containing 136 total intensity images. 
Analysis of the image sequence reveals a sharp jet-position-angle swing by > 100° during [2004,2006], as 
viewed in the plane of the sky, that we interpret as the crossing of the jet from one side of the line of sight 
to the other during a softer and longer term swing of the inner jet. Modulating such long term swing, our 
images also show for the first time a prominent erratic wobbling behavior of the innermost ~ 0.4 mas of the 
jet with fluctuations in position angle of up to ~ 40° over time scales ~ 2 yr. This is accompanied by highly 
superluminal motions along non-radial trajectories, which reflect the remarkable non-ballistic nature of the jet 
plasma on these scales. The erratic nature and short time scales of the observed behavior rules out scenarios 
such as binary black hole systems, accretion disk precession, and interaction with the ambient medium as 
possible origins of the phenomenon on the scales probed by our observations, although such processes may 
cause longer-term modulation of the jet direction. We propose that variable asymmetric injection of the jet 
flow; perhaps related to turbulence in the accretion disk; coupled with hydrodynamic instabilities, leads to the 
non-ballistic dynamics that cause the observed non-periodic changes in the direction of the inner jet. 

Subject headings: Galaxies: active - galaxies: jets - BL Lacertae objects: general - BL Lacertae objects: 
individual (OJ287) - radio continuum: galaxies - polarization 



1. INTRODUCTION 

The most luminous long-lived sources of radiation in the 
cosmos — active galactic nuclei (AGN) — are powered by gas 
falling from an accretion disk onto a super-massive black hole 
(BH, > 10 7 M Q ) at their center. AGN classified as blazars are 
characterized by wild variability of flux of non-thermal radi- 
ation from radio to 7-ray frequencies. Members of this class 
include BL Lacertae objects (BL Lacs) and, at higher emitted 
power, flat spectrum radio quasars (FSRQ). The remarkable 
properties of blazars include superlumin al apparent motions 
that can exceed 40 c dJorstad et al.ll2005l) . as well as substan- 
tial changes in flux and linear polarization on time scales as 
short as minutes. These phenomena are thought to be caused 
by relativistic jets of highly energized, magnetized plasma 
that are prop elled along the rotational p oles of the BH-disk 
system (e.g., Blandford & Znaiek 1977). AGN jets pointing 
within < 10° of our line of sight beam their radiation and 
shorten the variability time scales to give blazars their extreme 
properties. 

High-resolution long-term very-long baseline inter- 
ferometric (VLBI) monitoring of blazars is revealing a 
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growing number of cases where the observed (in pro- 
jection) jets undergo rapid swings of position angle in 
their innermost regions while remai ning collimated out 
to kiloparsec scal e s or farther (e.g.. iHoman et al] I2003t 
I Stirling et all 120031; rTatevama & KingharrJ l2004t iBach et all 
20051; lLobanoy & Roland! 120051; IMutel & Dennl 120051; 
ISavolainen et alJ[2006tlAgudo et al.ll2007h . Time scales rang- 
ing between ~ 2 and > 20 years and projected oscillations 
of the jet position angle with amplitudes fr om ~ 25° to 
> 100° are typical for the reported cases (see lAgudoll2009l 
for examples). This phenomenon — which has previously 
been termed jet precession, although we prefer wobbling 
given the doub ts about its periodicity in some sources (e.g., 
IMutel & Dennl 120051) — has been reported so far only in 
the innermost (parsec-scale) regions of jets. This already 
implies that the origin of jet wobbling might be a mechanism 
intimately tied to the origin of the jet. Bends and helical 
structures in AGN are also claimed to be produced by changes 
in direction near the base of the jet (e. g .|Dhawan et al.ll 1998h 
iLister et alJ 120031; lHardeeet alJl2005l lAgudo et all 12006a). 
which may be related to the wobbling phenomenon as well. 

During the past few decades, the BL Lac object OJ287 
(z = 0.306) has shown quasi-periodic double-peaked opti- 
cal flares every ~ 1 2yr, most recently in [2005,2008] (e.g., 
iVillforth et aDl20loh . This quasi-periodicity has been used to 
support binary BH models, as well as other scenarios such 
as accretion disk or jet instabilities (see IVillforth et alJl2010l 
for a review of these d ifferent scenarios). Consisten t with the 
binary BH hypothesis, Tatevama & Kinghami (120041) reported 
jet wobbling in a 3.5 cm long-term Very Long Baseline Array 
(VLBA) study from 1994 to 2002. They interpreted this as 
the result of ballistic precessio n of the jet, also with a peri- 
odicity of ~ 12 yr. Ba sed on thelValtonen et all ([2008) binary 
BH model for OJ287. |yaitonen et al.l d201 ll) have proposed a 



2 



Agudo et al. 



new model for wobbling of the jet in this source with a peri- 
odicity of ~ 120yr and a modulation of the jet position angle 
ona~ 12 yr time scale driven by changes in orientation of the 
inner accretion disk of the proposed primary black hole. 

To re-evaluate the origin of jet wobbling in general, and 
in OJ287 in particular, we present here the results of a de- 
tailed 7 mm VLBA study of OJ287 from 1995 to 2011 that 
includes 136 VLBI images. Owing to the reduced effect of 
synchrotron self-absorption in blazar jets at these short wave- 
lengths, compared with centimeter-wave radio bands, and the 
ultra-high (~ 0.15 milliarcsecond) angular resolution of the 
images, we scrutinize the inner jet regions in much finer detail 
than is possible at centimeter wavelengths. We adopt hereafter 
the standard ACDM cosmology with Ho=71 km s" 1 Mpc -1 , 
VIm = 0.27, and Vl\ = 0.73. Under this assumption, at the 
redshift of OJ287, 1 milliarcsecond (mas) translates to a pro- 
jected distance of 4.48 pc, and a proper motion of 1 mas yr" 1 
corresponds to a superluminal speed of 19.09c. 

2. OBSERVATIONS, DATA REDUCTION, AND 
MODELING 

We present 136 total-intensity VLBA images of OJ287 at 
7 mm obtained from January 1995 to June 2011; see Ta- 
ble Q] and Fig. Q] Linear polarization 7 mm VLBA im- 
ages from Dec ember 2010 to Jun e 2011 are also shown 
in Fig. [2] (see Ag udo et al.l 1201 laL for polarization maps 
before Dec. 2010). Within the time range 1997.86- 
2001.95, most of the data were obtained on a monthly 
basis under a program that mon itored the compact struc- 
ture of the r adio galaxy 3C 120 dGomez etalj|2000l |200ll 
120081 120111) . where OJ287 was one of the main calibra- 
tors. During the time frame 2006.38-2011.45, most of 
the observations were performed on a roughly monthly 
basis under the Boston University 7 mm blazar monitor- 
mg program JiyArcangelo et al. 2009; Larion ov et all 12 009; 
Marscheret al. 2010 ; [Jorstad et alJl2010HAgudo et alj|201 lal 
Ichatterjee et al] 1201 lb lAgudo et alJ 1201 lbl) . The data from 
these programs are also supplemented by OJ287 observations 
from several other in dependent 7 mm VLBA programs dat- 
ing back to 1995 (e.g.Uorstad et al.ll200H lAgudo et al.ll 2006a: 
iPiner et al.ll2006t lAgudo et alJl2007t lAgudd l2009). To cover 
gaps in the time coverage of the source, we obtained data at 
eight observing epochs from the VLBA archive. 

The data were reduced following the standard method for 
high frequency VLBA polarimetric dat a describ ed in, e.g., 
Jorsta d et al] (120051) : lAgudo et al] (l2006al) . and iGomez et alJ 
(120111) . For further details on the data reduction procedure 
followed for each particular data set, we refer to the corre- 
sponding publications indicated in Table Q] 

To make a simpler representation of the brightness distribu- 
tion in our images, we have performed model fits in which 
the emission features are represented by circular Gaussian 
emission components, with the number of components kept as 
small as possible. To estimate the uncertainties o f the model 
fit par ameters, we follow the prescription of Jorst ad et alJ 
(2005). Our fits typically include no more than two or three 
compact components within the innermost ~ 0.2-0.4 mas of 
the jet (Tabled. We use the position of these components in 
the inner ~ 0.4 mas of the observed jet to define the jet posi- 
tion angle (JPA) in OJ287 as observed from Earth (see Fig. [3). 

3. RESULTS AND ANALYSIS 
3.1. The Innermost VLBI Jet Feature 



We define the "core" of the brightness distribution of OJ287 
as the innermost jet feature visible in our 7 mm VLBA im- 
ages. Given tha t the core is not us ually associated with mov- 
ing blobs (e.g. lJorstad et alJl2005l) we assume that it is a sta- 
tionary feature in the jet of OJ287. Identification of the core in 
a multi-epoch kinematic VLBI study like ours — where self- 
calibrated VLBI images do not preserve the absolute posi- 
tional information — is important to define a kinematic refer- 
ence point of the system. 

For the time period from the beginning of 1995 to mid- 
2004, the identification of the core is clear. The bright jet 
feature at the eastern end (northeastern end at some epochs) 
of the brightness distribution, i.e., model-fit component C (see 
Fig- UK from which supe rlum inal features emerge and travel 
downstream (see Section 1331 ), corresponds to the core. How- 
ever, the identification of the core is more complicated dur- 
ing the time range from mid-2004 to mid-2011. Whereas 
the out er knots (beyond 0. 5 mas from component C, labelled 
CO in lAgudo et alj|201 lah continue propagating toward the 
southwest, the region within the innermost — 0.4 mas behaves 
in a completely different way than observed before. The mo- 
tion of knots emerging from the core suddenly shifted to the 
north (and later to the northwest) direction after mid-2004. 
The new direction is nearly perpendicular to the typicaljet 
position angle during the 1995-2003 time range (see Fig . |Tj> . 

The identification of the core with component C is justified 
by: (i) the detection of superluminal motion away from com- 
ponent C, toward jet feature a, and then outward (see com- 
ponent nfj; (ii) the progressive decrease of jet intensity to- 
ward th e west of component a (labelled CI in Agu do et alJ 
d201 lal) . see Fig. []} at epochs after mid-2010; (iii) the con- 
firmation of a continuous jet flow tracing the arc C-a-n af- 
ter mid-2010 in both total and polarized intensity and in po- 
larization angle (which remains parallel to the local jet axis; 
see Fig. 12); and (iv) the smaller mean size and higher mean 
observed brightness-temperatur^B of component C from fall 
2004 (0(C) = (0.02 ±0.02) mas and 7b(C) - 2.7 x 10 12 K, re- 
spectively.) as compared to that measured for component a 
(0(a) = (0.05 ±0.03) mas and r B (a) - 1.1 x 10 12 K) on the 
same dates. All these lines of evidence together suggest that 
the jet flows downstream from component C toward compo- 
nent a, and at epochs after mid-2010 toward components m 
and n. This unambiguously identifies component C with the 
innermost visible jet feature in our 7 mm images, i.e., with the 
"core" of OJ287 as defined above. 

3.2. Jet Wobbling in OJ287 

Figures [T] and [3] show clear clockwise swing of the jet in 
OJ287 over the time range between the beginning of 1995 
(with JPA- -100°) and mid-2004 (with JPA- -150°). This 
trend is modulated by shorter term jet wobbling on time scales 
— 2 yr (as measured from local JPA maxima and minima) and 
relatively small amplitude, > 25°. During the [1995,2004] 
time range, this phenomenology is not necessarily inconsis- 
ten t with the interpretatio n of str ictly periodic jet precession 
by Tatevama & Kingham (2004), as long as the short-term 
wobbling on fine angular scales seen in our data is superposed 

8 this is also clearly seen in the evolution of the total inten- 
sity of the jet in an animation at URL https://w3.iaa.es/ 
^iagudo/ research/O J 2 87 /O J 2 87_2 1 0_e jection . mov, 

as well as in Figs.[f]and|2]and TablefH 

9 Th e o bserved brightness tem perature was computed as in Marsc her et alJ 
1 1 979 ) and Agudo et al. 1 2006b ) 
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FIG. 1 . — Sequence of 7 mm VLBA maps of OJ287 from 1995 to 201 1, convolved with a circular Gaussian beam with FWHM = 0.15 mas (represented by the 
circle in the left-bottom corner). Contours symbolize the observed total intensity with levels in factors of 2 from 0.1 to 51.2 % plus 90.0 % of peak= 9.018 Jy/beam. 
A movie corresponding to this sequence of images can be found at URL https://w3.iaa. es/~iagudo/research/OJ2 8 7 /OJ2 8 7LT . mo v. 



on long-term, sinusoidal modulation (with period ~ 12 yr) as 
inferred by these authors on coarser angular scales. In addi- 
tion to the wobbling of the jet, our images feature superlumi- 
nal knots that propagate down the jet at speeds ranging from 
~ 2c to ~ 13c (see Tabled. 

In mid-2004, the 7 mm brightness distribution of OJ287 
is dominated by emission from the nearly point-like core 
(Fig-ID- A few months later, the core region starts to expand 
essentially in the north-south direction. By the end of 2005, 
and even more so by the end of 2007, this expansion allows us 
to discern two separate emission regions oriented in the north- 
south or northwest-southeast direction. This relatively sudden 
change in the structure of the jet leads to an extremely sharp 
swing of JPA in mid-2004, as shown in Fig. [3] Given that this 
sudden jump in JPA is produced during a relatively short time 
period (< 1-2 yr) when the dominant emission structure is ex- 
tremely compact, we propose that the event is produced by a 
progressive change of orientation with respect to the observer. 
The innermost jet axis should have passed from one side of the 
line of sight (before mid-2004), to essentially along the line of 
sight (in mid-2004), and finally to the other side of the line of 
sight (after mid-2004). This can explain the sharp, relatively 



short change of JPA in mid-2004, as projected on the sky in 
our images (Figs.[T]and[3]l. This interpretation is further sup- 
ported by the sharp peak of core emission at about the time of 
the sharp JPA swing in mid-2004 (Fig. [3]) owing to enhanced 
Doppler boosting of the jet emission when the it points closer 
to the obser ver's line of sight. 

Following lJorstad et alJ (120051) . we have estimated the vari- 
ability Doppler factors (<5 var ), viewing angles (# var ), and 
Lorentz factors (r var ) from the decaying flux density patterns 
and superluminal speeds of well-defined and bright compo- 
nents in our images of OJ287 (TableHJ). To perform estimates 
of these physical parameters that are as reliable as possible, 
we select only jet features with well-determined positions and 
maximum apparent speeds. These are associated with iso- 
lated jet features or with a leading perturbation, rather than to 
trailing features, which are known to display intrinsic pattern 
speeds mu ch lower than tho se of the corresponding leading 
features ( Agud oet al.ll200l1) . We also require the model-fit 
components to be identified over more than five epochs, with 
peak (mean) flux densities higher than 200 mJy (90 mJy), and 
with clearly decaying flux density patterns. The scatter in 
our time-dependent estimates of viewing angle A6\, al - < 2.7, 
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FIG. 1. — Continued. 

which we primarily attribute to small uncertainties in the de- 
termination of 6> var , does not allow us to verify the proposed 
hypothesis of the crossing of the jet from one side of the line 
of sight to the other. However, the small values of the viewing 
angle, 6» var <E [0.7,3.4], during the 1997-2010 time range fa- 
vors this hypothesis, which is more likely for jets with close 
to zero. 

The new brightness distribution of OJ287 after the extreme 
JPA change in mid 2004 (see Fig. [3]) also shows clear wob- 
bling of the jet from mid-2004 to mid-201 1. Such long term 
wobbling, which is essentially tracked by the JPA defined by 
the core and the nearby, bright jet feature a, causes the jet 
to rotate in the plane of the sky from ~ 0° at the beginning 
of 2005 to ~ -50° at the beginning of 2011. Superposed 
on this long-term JPA variation after mid-2004, we also find 
clear short-term JPA variability. Although the amplitude of 
these short-term fluctuations in the 2004.5-201 1.5 time range 
(~ 40°) seems to be larger than that from 1995 to 2004.5, the 
time scale is roughly similar in both cases (~ 2 yr). 

3.3. Non-Ballistic Superluminal Motion 

We use the model fits performed for every image to identify 
all jet features that persist over several epochs or more (Fig.Q] 




and Tables|2]and[3]). From these identifications, we derive the 
trajectories of every prominent jet feature (with average flux 
density < S >> 40mJy). Among all identified components at 
several epochs we find both stationary and moving features, 
the latter showing both rectilinear (or slightly bent) and non- 
radial superluminal motion (see Table [3]). For rectilinear tra- 
jectories (with superluminal speeds in the range 2-13 c), we 
calculate linear fits of r (distance to the core) versus t (time) 
to estimate average apparent speeds and ejection times (to). 
For non-radial superluminal trajectories (showing apparent 
speeds up to ~ 15 c) we follow the second-order polynomial 
method of Agu do et"aTI (120071) to fit the trajectories and de- 
termine the kinematic properties of the corresponding jet fea- 
tures. Although superluminal non-r adial motions hav e been 
reported before, e.g., in NRAO 150 dAgudo et alj|2007l) . this 
is the first time that superluminal non-radial motions > 10 c 
have been reported in a blazar (jet feature 1, TableO. 

We find highly non-radial superluminal motions in the com- 
pact jet of OJ287 only after mid-2004, when the projected jet 
axis appears to be sharply bent. The moving jet features dur- 
ing the 2004.5-2011.5 time range (i.e., components a, I, m, 
and «, Table |3]l follow highly bent trajectories. This demon- 
strates the non-ballistic motion of the jet fluid on the scales 
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FIG. 1. — Continued. 

covered by our VLBA observations during this time range. 
The fact that only rectilinear projected trajectories are ob- 
served during the time range when the jet exhibits a more 
rectilinear shape, i.e., 1995-2004.5, does not rule out that even 
during this time range the fluid is non-ballistic in three dimen- 
sions. 

4. DISCUSSION 

4. 1 . Erratic versus Periodic Jet Wobbling 

Neither the short-term nor long-term variations in JPA dis- 
played by Fig.[3]correspond to a simple pattern that we can re- 
late to regular or periodic motion, at least over the time scales 
covered by our observations (~ 16 yr). Instead, they seem to 
be related to some erratic process acting on short observed 
time scales < 2 yr with amplitudes < 40° (see JPA swings 
during 1995, 2007, 2009, and 201 1 in Fig.©. 

Although based on poorer time and spatial sampling than 
that employed in our study, slower JPA changes than those 




reported here have been suggested before on larger pro- 
jected scales (tens of parsecs) than those probed by our 
7 mm o bservations (a few p arsecs); see lTatevama & Kingham 
(2004); Mooretal. (2011), who used 3.5 cm observations, 
and IValtonen et alj (1201 ll) . who employed VLBI observa- 
tions at 2 to 6 cm. VLBA observations at similar ra- 
dio wavelengths ( 2 cm) compiled by the MOJAVE teairFI 
dLister et al.ll2009[) over more than 10 y e ars also supports th e 
results of iTatevama & Kingham! ( 120041) : iMoor etalJ d2011h : 
Valton en et al.l d2011[). However, the periodic precession 
model proposed by|Tatevama & Kingham (2004) predicted a 
jet structural position angle ~ -100° in 2007, whereas we 
measure a drastically different inner JPA of ^ 0°. Hence, 
the jet structure shown in Fig. [T] does not follow Tateyama 
& Kingham's 12 yr periodic precession model. The JPA ob- 

10 www .physics . pur due . edu/astro/MO JAVE/ animated/ 
0851+202 . i .mpg 
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FIG. 1. — Continued. 



served by Mo6r et alJ (f20lTh is also > 100° off our 7 mm mea- 
surements on similar dates. Moreover, the JPA evolution on 
the scales observed through our observations (Fig. [3]) conflicts 
with t he JPA evolution model presented by Valton en et al.l 
( 1201 11) . This model predicts a JPA value in the range -1 10° to 
-150° at the beginning of 2009, whereas we measure ~ -10°. 

This apparent inconsistency could be reflecting a (perhaps 
periodic) long-term JPA variation in the jet on larger scales 
than those probed by our 7 mm VLBA observations, together 
with more erratic behavior modulating the long-term JPA 
variations in the innermost jet regions probed by our millime- 
ter monitoring. It thus appears that a single phenomenon can- 
not explain both the JPA variations observed on the scales of 
tens of parsecs, and the much shorter term (~ 2 yr) JPA varia- 
tions reported in this work. This implies that processes other 
than those producing the long-term JPA variations at large jet 
scales need to take place near the jet origin to modulate the jet 
direction on shorter temporal and spatial scales. 
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4.2. Possible Origin of Observed Jet Wobbling 

Different explanations have been proposed for wobbling 
in blazar jets. Orbital motion of compact objects in a bi- 
nary black hole system has be en suggested as a poss ibility 
to explain bends in AGN jets (Beg elman et al] 11980). Al- 
ternatively, accretion disk precession could arise from ei- 
ther a companion super-massive black hole or another mas- 
sive object inducing torques in the accretion disk of the pri- 
mary (e.g.JLister et a l. 2003), or by the Lense-Thirring effect 
(Bardeen & Pettersonlll975l) in the presence of a Kerr black 
hole with its rotational axis misaligned with regard to that of 
the accretion disk (e.g. Caproni et al. 2004). However, these 
processes cannot, by themselves, explain the short time scales 
(~ 2 yr) and erratic nature of the JPA changes reported here. 

Pressure gradients in the ambient external medium are also 
not, by themselves, a feasible explanation for the short-term 
JPA variations reported here. They might, however, explain 
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FIG. 3. — Top: Time evolution of the 7mm jet position angle (JPA) of 
OJ287. Different colors indicate the various jet features whose positions are 
used to estimate the JPA; see also Fig [T] Bottom: 7 mm total flux density 
evolution of the core. 

JPA variations that modulate the behavior of the innermost jet 
of OJ287. We discuss this possibility in the following subsec- 
tions. 



0.0 -0.5 -1.0 -1.5 
Relative Right Ascension (mas) 



FIG. 2. — Similar to Fig.[T]but restricted to the period December 2010 to 
June 2011. In addition, linearly polarized intensity (color scale) and polar- 
ization electric vector position angle distributions (short black sticks) are su- 
perimposed at each observing epoch. 

the smoothly be nt jet structure observed o n scales of hundreds 
of kiloparsecs (Marscher & Jorstad 201 1). 

A geometrical model based on non-ballistic motions of 
blobs and misidentification of components has also been pro- 
posed as_a_possible explanation of JPA variations in general 
(Gong 2008). However, we can exclude misidentification in 
our observations of OJ287 owing to the good time sampling 
of our monitoring. 

Other kinds of more erratic disk or jet instabilities 
(e.g., similar to those prod uced by turbulent accretion; 
McKinney & Blandford 2009) have not been extensively ex- 
plored previously as an explanation of jet wobbling in blazar 
jets. Such instabilities could explain the apparently erratic 



4.3. Jet Wobbling Produced by Variations of Flow Injection 

Since the fast and erratic JPA modulations that we re- 
port in this paper are observed in the innermost jet regions, 
we propose that time variations occur in the injection of 
the jet flow at the origin of the jet, perhaps resulting from 
changes in the particle density or magnetic field configura- 
tion caused by turbulence or irregularities in the accretion 
process. Such variations, combined with the axial rotation 
of the jet material in the presence of significant toroidal mag- 
netic field components, can produc e non periodic JPA varia- 
tions dMcKinnev & B landford 2009). If transverse gradients 
in pressure and velocity at the jet injection persist for extended 
time intervals, th ey should generate bent jet structures (e.g., 
lAlov et aT1l2003l) . that — amplified by projection effects from 
small mean jet viewing angles — may produce sharp apparent 
bends such as those observed in OJ287 (see Fig. [2]). 

4.4. Large Apparent JPA variations and Sharp Jet Bends 

The small viewing angle of the jet axis in some sources 
allows us to see the jet from inside its cone, in which case 
asymmetric changes in the jet can cause it to appear to cross 
the line of sight. The sharp jump in JPA in mid-2004 from 
JPA w -150° to JPA « 0° seen in o ur 7 mm sequence of im- 
ages supports this scenario (Section [372]). This requires only 
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minor changes of the jet structure provided that the mean 
viewing angle of the jet is close to 0°, as we show for OJ287 
(see Table|4|i. 

Projection effects from very small viewing angles can also 
create the illusion of very sharp bends from slightly bent 
jet structures, as is the case of the main bend in the jet of 
OJ287 from north-northeast to southea st (see Fig. [2). This 
bend could result from jet instabilities ( Hard eel20 1 UlPeruchol 
causing the jet flow to curve. Such instabilities can 
be triggered by v ariat ions at the jet injection as those men- 
tioned in Section 14.31 Moreover, it is reasonable to expect 
pressure gradients in the external medium to play an impor- 
tant role in re-directing the northeasterly flow toward the pre- 
vious JPA of ~ -1 10° at the edge of the channel carved out 
by the jet over thousands of years. Such bending would be ac- 
complished via multiple oblique shocks at the outer boundary 
of the jet. This in turn would be accompanied by compression 
of the magnetic field and acceleration of relativistic electrons, 
which would (1) enhance the synchrotron emission near the 
outer edge and (2) align the magnetic field near the outer edge 
so that it follows the curvature of the jet. The polarization map 
of 2011-04-21 (and possibly 2011-06-12; see Fig.© conform 
with these expectations, although the angular resolution in the 
north-south direction is inadequate to confirm the predicted 
edge-brightening. 

4.5. Relation to Long Scales 

If inhomogeneous jet injection occurs randomly, i.e., 
with no systematic temporal pattern, its effect on the jet 
flow in terms of instabilities would decrease with dis- 
tance and within the jet cross section. Indeed, the mean 
long-term position angle o f the compact jet found by 
Tatev ama & Kinghaml J20 04) is ~ -107°, very similar to the 
value of -109° of the kiloparsec -scale jet within 8" of the nu- 
cleus dMarscher & Jo rstad 201 1). This implies that only short 
wavelengths are excited, perhaps in the form of high-order 
body modes. Such short wavelengths grow rapidly in am- 
plitude but are inefficient in jet mixing and disruption in the 
non-linear phase, in contrast with instabilities that generate 
large deformati ons on the jet surface, i.e., first-order b o dy and 
surface modes dXu et alJlioOOt iPerucho et alJ l2004alH 120051 
120071 120 lOt lHardee & Eilekll201 11) . Thus, the small-scale na- 
ture of the observed changes in our 7 mm VLBA data do not 
imply disruption of the jet on larger scales despite the high 
apparent amplitudes, which are accentuated by projection ef- 
fects. 

5. CONCLUSIONS 

Our 7 mm VLBA study has revealed a sharp jet-position- 
angle swing by > 100° during [2004,2006], as viewed in the 
plane of the sky, that we interpret as the crossing of the jet 
from one side of the line of sight to the other during a softer 
and longer term swing of the inner jet. Modulating such 
long term swing, our images also show that the innermost 
~ 0.4 mas region of the 7 mm jet in OJ287 wobbles follow- 
ing an erratic pattern with amplitudes up to ~ 40° on time 
scales of ~ 2 yr. The latter phenomenon, reported here for 
the first time, rules out scenarios such as the putative orbit 
of the supermassive black hole responsible for the produc- 
tion of the jet in a binary system, accretion disk precession 
induced by either a companion compact object or the Lense- 
Thirring effect, and jet interaction with the external medium 
as possible origins of the phenomenon. This is based on the 



long time scales related to these scenarios compared with the 
~ 2 yr time scale and to the apparently erratic nature of the 
observed wobbling in OJ287. The dense time sampling of our 
observations also allows us to rule out misidentification of jet 
features across epochs as a possible explanation. This leads 
us to propose a new possibility to explain the non-periodic jet 
wobbling reported here through asymmetric fluctuations in in- 
jection of plasma flow into the jet coupled with instabilities. 
This could be related to variations in the accretion process and 
in the disk-jet connection. 

Both our observations an d those from other VLB I programs 
at cen timeter wavelengths dMoor et alJ 1201 lb IValtonen et alJ 
1201 1[) rule out the ~ 12 yr periodi c precession model proposed 
by iTatevama & Kingham ( 120041) . which predicted different 
evolution of th e JPA than observed. I n contrast, the 120 yr pe- 
riodicity of the Valto nen et alJ (1201 ll) precession model seems 
able to reproduce the JPA evolution observed at centimeter 
wavelengths, although the 12yr modulation that their model 
imparts on the predicted JPA changes does not clearly fol- 
low their observed JPA variations during the [2005,201 1] time 
range. Such 12 yr modulation mig ht be better rep roduced by 
the 3.5 cm JPA observations by Moor et al. (201 1), but more 
data is still needed to clarify the consistency between new Val- 
tonen' s JPA evolution-model and the centimeter VLBI jet be- 
havior. 

Our interpretation of the apparent mismatch between the 
JPA evolution reported in this paper (i.e., on scales of a few 
parsecs) and those observed at centimeter wavelengths (tens 
of parsecs) leaves open the possibility that different instabil- 
ity modes operate on these two spatial scales. This allows for 
the JPA on the inner scales to behave in an erratic way, while 
on the larger scales it can be driven by instabilities coupled 
to modes with periodic behavior. Periodic models based on 
binary black hole systems are thus not ruled out by our ob- 
servations, unless future centimeter VLBI observations of the 
jet in OJ287 show more erratic JPA variations similar to those 
seen at millimeter wavelengths. 

The detailed study of the evolution of the JPA of OJ287 pre- 
sented here, and the main conclusions drawn above, were only 
possible with the ultra-high resolution provided by the VLBA 
at 7 mm and long-term monitoring of OJ287 with dense time 
sampling. Studies similar to that presented here can be car- 
ried out for a number of blazars that have been monitored in- 
tensively with millimeter VLBI during the past decade. Per- 
forming these studies will be crucial to evaluate how common 
erratic wobbling is in the jets of AGN, and to draw general 
conclusions about the possible relation of this phenomenon to 
the accretion process. 



This research has made use of the M OJAVE database main- 
tained by the MOJAVE team dLister et alJl2009h . The VLBA 
is an instrument of the National Radio Astronomy Observa- 
tory, a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
The research at the IAA-CSIC is supported in part by the 
Ministerio de Ciencia e Innovacion of Spain, and by the re- 
gional government of Andalucfa through grants AYA2010- 
14844 and P09-FQM-4784, respectively. The research at 
Boston University was funded by US National Science Foun- 
dation grant AST-0907893, NASA grants NNX08AJ64G, 
NNX08AU02G, NNX08AV61G, and NNX08AV65G, and 
NRAO award GSSP07-0009. 
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TABLE 1 

Observing log and image information. (Truncated table). 



Epoch 


tint 


Bit 






^peak 


Noise 


Ref. 




min 


samp 


MHz 


Jy 


Jy/beam 


mJy/beam 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 



1995-01-05 


68 


1 


32 


1.01 


0.87 


2.5 


1 


1995-04-18 


25 


1 


32 


1.29 


0.93 


3.1 


1 


1995-06-04 


34 


1 


32 


1.15 


0.93 


1.8 


1 


1995-10-15 


14 


1 


32 


2.57 


1.93 


4.4 


1 


1996-05-03 


30 


1 


64 


1.18 


0.90 


1.6 


1 


1996-1 1-24 


35 


1 


64 


1.19 


0.85 


1.3 


1 


1997-07-30 


35 


1 


64 


1.45 


1.09 


0.8 


1 


1997-11-10 


30 


1 


64 


1.48 


1.06 


0.6 


2 


1997-12-11 


32 


1 


64 


1.26 


0.81 


0.6 


2 


1998-01-11 


32 


1 


64 


1.58 


1.13 


1.3 


2 


1998-02-07 


32 


1 


64 


1.83 


1.38 


0.7 


2 



REFE RENC ES. — (1) IJorstad et alj 120011): ( 2) Gomez et al. 1200 1)); (3) | Agudo et alj 
:2006a); (4) Bach et al. i 2006); (5) Gomez et al. i2008); (6) IGomez : et al) 1201 II) ; 

(7) iPiner et alj (20061): (8) da ta from VLB A arc hive; (9) lAgudoetalJ 120071) : (10) 
iD'Arcangelo et alj 120091) : (1 D IAgudo etai] 120101) : q2) IAgudo et alj 1201 lal) : (13) this 
paper. 

NOTE. — Columns are as follows: (1) observing epoch, (2) total integration time, (3) 
number of bits used for signal digitalization sampling, (4) observing frequency bandwidth, 
(5) total integrated flux density, (6) peak flux density, (7) noise level of the resulting image, 

(8) reference where data was published. 
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TABLE 2 

Model fit parameters. (Truncated table). 



Comp 


S 


r 


6 


FWHM 




mjy 


mas 


O 


mas 


(1) 


(2) 


(3) 


(4) 


(5) 




1995-01-05 


C 


902 ±60 


0.00 ± 0.00 


0±0 


0.022 ±0.001 


A 


96±7 


0.23 ±0.03 


-83 ±7 


0.176 ±0.009 


B 


55±4 


0.62 ±0.03 


-88 ±2 


0.243 ±0.012 


1995-04-18 


C 


1034 ±69 


0.00 ± 0.00 


0±0 


0.061 ±0.001 


A 


210±15 


0.24 ±0.03 


-95 ±7 


0.130 ±0.006 


B 


54±3 


0.78 ±0.03 


-94 ±2 


0.259 ±0.013 


D 


63±4 


1.09 ±0.03 


-90 ±1 


0.186 ±0.009 


E 


91 ±10 


1.40 ±0.25 - 


-92 ±10 


0.5 10 ±0.051 


1995-06-04 


C 


1034 ±69 


0.00 ±0.00 


0±0 


0.049 ±0.001 


A 


41±4 


0.30 ±0.01 - 


-118±1 


0.000 ±0.001 


D 


142 ±17 


1.11 ±0.29 - 


-91 ± 14 


0.570 ±0.057 



NOTE. — Columns are as follows: (1) jet feature, (2) flux density, (3) 
distance to the core, (4) position angle, (5) FWHM. 
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TABLE 3 

Properties of model fit components. 





Comp 


<S> 


<e> 


< r > < Vapp > 


< V ;lp p > 


to 


^ .,nonradial \ 
^ "app ■* 


^ ^ipp > 


<0> 






mJy 


O 


mas mas/yr 




yr 


mas/yr 




°/y 




(1) 


(2) 


(3) 


(4) (5) 


(6) 


(7) 


(8) 


(9) 


(10) 












Quasi-stationary or slowly moving 


features 








A a 




294 ± 182 


-1 13 ± 8 


0.28 ±0.06 












F 




47 ±25 


-98 ±6 


1.38±0.11 












G b 




77 ±41 


-110±2 


1.02±0.06 












L c 




174 ±148 


-102 ±15 


0.15±0.04 












e 




158 ±111 


-89 ±19 


0.26 ±0.07 












f 




400 ±228 


-70 ±16 


0.19±0.03 












g 




89 ±36 


-127 ±9 


0.74 ±0.03 












k 




200 ±138 


-57 ±12 


0.18±0.02 












Superluminal rectilinear features 


H d 




93 ±63 


-116±6 


0.43 ±0.08 0.32 ±0.08 


6.0 ±1.5 


1997.71 ±0.25 








J e 




119±99 


—111 ±3 


0.57 ±0.16 0.40 ±0.07 


7.7 ± 1.3 


1998. 02 ±0.24 








K< 




37 ±24 


-113±8 


0.79 ±0.34 0.41 ±0.04 


7.7 ±0.7 


1999.22±0.18 








0« 




243 ±135 


-116±4 


1.19±0.26 0.50±0.01 


9.5±0.1 


2000. 16 ±0.03 








N 




288 ±152 


-122 ±2 


0.58 ±0.09 0.39 ±0.14 


7.5 ±2.6 


2000.76 ±0.33 








R 




91 ±16 


-123 ±1 


0.85±0.18 0.68±0.06 


13.0±1.2 


2001.92±0.12 








T 




97 ±46 


-124 ±6 


0.97 ±0.30 0.40 ±0.08 


7.6±1.5 


2003.22 ±0.38 








V 




221 ±314 


-130± 10 


0.51 ±0.44 0.54 ±0.05 


10.2±0.9 


2004.30 ±0.10 








b 




72 ±36 


-120±3 


1.17 ±0.13 0.46±0.15 


8.8±2.8 


2006.34 ±0.48 








d 




71 ±46 


-126 ± 14 


0.60 ±0.24 0.24 ±0.01 


4.5 ±0.2 


2004.77 ±0.13 








h 




120 ±50 


-98 ± 22 


0.25 ±0.09 0.30 ±0.02 


5.7 ±0.4 


2006.78 ±0.06 








j 




157 ±49 


-126 ±7 


0.81 ±0.39 0.35 ±0.02 


6.7 ±0.4 


2007.85 ±0.12 








Superluminal bent features 


a 




2293 ±1631 


-22 ±23 


0.17 ±0.06 0.07±0.01 h 


1.4±0.1 h 




0.07±0.01 h 


1.4±0.1 h 


-19±l h 


1 




444 ± 404 


-62 ±21 


0.33 ± 0.07 0.94 ±0.21 


17.9±4.0 




0.78 ±0.29 


14.9 ±5.5 


-145 ±40 


m 




463 ±284 


-75 ±9 


0.35 ±0.04 0.62 ±0.14 


11.9±2.6 




0.52±0.18 


10.0±3.5 


-87 ±21 


n 




531 ±149 


-67 ± 13 


0.40 ±0.1 3 0.54 ±0.02 


10.3 ±0.5 




0.22 ±0.09 


4.3 ±1.7 


-40±4 



NOTE. — Columns are as follows: (1) jet feature, (2) average flux density, (3) average position angle, (4) average distance to the core, (5) average proper motion in mas/yr, (6) 
average proper motion in units of c, (7) estimated ejection time if any, (8) average proper motion in non radial direction in mas/yr, (9) average proper motion in non radial direction 
in units of c, (10) average angular speed. 

a This jet feature corresponds to stationary component A2 as identified by Jorstad et al. 1 2005), and CI as identified by Jorstad et al. ( 2001). 
k This jet feature corresponds to stationary component A3 as identified by Jorstad et al. i 2 0051) . 
c This jet feature corresponds to stationary component Al as identified by Jorstad et al. 1 2 0051) . 

This feature may correspond in some epochs to superluminal knot B4 identified by Jorstad et al. 1 2005). We prefer our identification owing to our better time sampling during 
1998 though. 

e This feature corresponds to superluminal knot B5 identified by Jorstad et al. 1 2005), although they reported faster superluminal speed (12 ± 2 c) perhaps owing to their better time 
sampling during 1999, hence producing better kinematic estimates. 

^ This feature may correspond to superluminal knot B6 identified by Jorstad et al. 1 2005). We prefer our identification owing to our better time sampling during 2001 though. 
g This feature corresponds to superluminal knot B7 reported by Jorstad et al. i 2005 ). 
Computed only for the time range from 2008-10-22 to 2010-10-24, when the bent trajectory of component a is better defined and reaches superluminal non-radial speeds. 



TABLE 4 

Physical parameters of well defined 
superluminal components 



Comp 


to 


5var 


Ovai 


l-Var 




yr 




deg. 




(1) 


(2) 


(3) 


(4) 


(5) 


H 


1997.71 ±0.25 


15.5 


2.7 


8.7 


J a 


1998.61 ±0.09 


16.1 


3.4 


12.3 


O 


2000. 16 ±0.03 


33.3 


0.9 


18.0 


V 


2004.30±0.10 


22.6 


1.9 


13.6 


h 


2006.78 ±0.06 


25.2 


1.0 


13.3 


n 


> 2010.64 


36.2 


0.7 


19.1 



il We show the parameters computed by Jo rstad et"afl 
( 2005), who made a better determination of the kinemat- 
ics of this jet feature (see Table[3j. 



